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Abstract 

Three- and four-jet final states have been measured in photoproduction at HERA 
using the ZEUS detector with an integrated luminosity of 121 pb _1 . The results 
are presented for jets with transverse energy > 6 GeV and pseudorapidity 
\7f et \ < 2.4, in the kinematic region given by the virtuality of the photon Q 2 < 
1 GeV 2 and the inelasticity 0.2 < y < 0.85 and in two mass regions defined as 
25 < M n j < 50 GeV and M n j > 50 GeV, where M n j is the invariant mass of the 
n-jet system. The four-jet photoproduction cross section has been measured for 
the first time and represents the highest-order process studied at HERA. Both the 
three- and four-jet cross sections have been compared with leading-logarithmic 
parton-shower Monte Carlo models, with and without multi-parton interactions. 
The three-jet cross sections have been compared to an 0(aa 2 s ) perturbative QCD 
calculation. 
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1 Introduction 



In photoproduction at HERA, a quasi-real photon, emitted by the incoming electrons, 
interacts with the proton. Hard photoproduction [1-10] may be categorised at leading 
order (LO) as being either direct, if the photon interacts as a point-like particle, or 
resolved, if it fluctuates into a partonic system, and subsequently transfers only a fraction 
of its momentum to the hard interaction. 

Photoproduction collisions at HERA can lead to final states with multiple jets. Multi-jet 
events are of particular interest as they are produced by processes that are manifestly 
beyond LO in the strong coupling constant, a s . Presently, predictions from perturbative 
quantum chromodynamics (pQCD) pertaining to multi-jet final states in photoproduction 
are only available up to 0(aa 2 ), where a is the fine structure constant. 

The hadron-like structure of the photon in resolved processes gives rise to the possibility 
of multi-parton interactions (MPIs) at HERA. In the MPI picture, more than one pair 
of partons takes part in the hard interaction. A schematic of an MPI event is shown 
in Fig. [TJ The secondary scatters generate additional hadronic energy flow in the event, 
the topology and magnitude of which is poorly understood theoretically. Potentially, this 
energy flow may lead to the formation of jets and so MPIs may constitute a source of 
multi-jets in the final state. Multi-parton interactions have been studied before at the 
Tevatron [11-14], in dijet photoproduction events at HERA [15,16] and are expected to 
be prevalent at the LHC [17,18]. 

Three-jet photoproduction events have been studied before at HERA in the M$j > 
50 GeV [19] region of three-jet invariant mass. Multi-jet final states have also been 
studied in deep inelastic scattering (DIS) at HERA [20-23]. In this paper, differential 
cross sections are presented for the three-jet photoproduction final state, in a wider M 3j - 
region and with over seven times the luminosity of the previous ZEUS publication [19]. 
The four-jet photoproduction cross sections are presented here for the first time. Also 
examined is the description of the data by two Monte Carlo (MC) models and whether 
or not it is improved by the introduction of simulated MPIs. In addition, an 0(aa 2 s ) 
prediction has been compared to the three-jet data. 



2 Experimental conditions 

The data were collected using the ZEUS detector during the 1996 to 2000 running periods. 
In 1996 and 1997, HERA collided positrons with protons, with energies of E e = 27.5 GeV 

1 From now on the word "electron" is used as a generic term for electrons and positrons, unless stated 
otherwise. 
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and E p = 820 GeV, respectively, corresponding to a beam centre-of-mass energy, ^/s, 
of 300 GeV. From 1998 onwards the proton beam energy was raised to E p = 920 GeV 
(y/s = 318 GeV). Furthermore, a subsample of these later data consists of electron-proton 
collisions. The total sample corresponds to an integrated luminosity of (121 ± 2) pb _1 , 
of which 82 pb _1 were collected at yfs = 318 GeV. A detailed description of the ZEUS 
detector can be found elsewhere [24,25]. A brief outline of the components that are of 
most relevance to this analysis is given below. 

The high-resolution uranium-scintillator calorimeter (CAL) [26] consists of three parts: 
the forward, the barrel and the rear calorimeters. Each part is subdivided transversely 
into towers and longitudinally into one electromagnetic and either one (in the rear) or two 
(in the barrel and forward) hadronic sections. The smallest subdivision of the calorime- 
ter is called a cell. The CAL relative energy resolutions, as measured under test-beam 
conditions, are 0.18/ VE for electrons and 0.35/\^E for hadrons, with E in GeV. 

Charged particles are measured in the central tracking detector (CTD) [27], which op- 
erates in a magnetic field of 1.43 T provided by a thin superconducting coil. The CTD 
covers the polar-anglj§ region 15° < 9 < 164°. The relative transverse-momentum resolu- 
tion for full-length tracks is a PT /p T = 0.0058p T ©0. 0065©0. 0014/pr, with p T in GeV. Both 
tracking and calorimetry were used to reconstruct the transverse energy and direction of 
jets as described in Section [31 

The luminosity was measured from the rate of the bremsstrahlung process ep —>■ ejp, 
where the photon is detected in a lead-scintillator calorimeter [28] placed in the HERA 
tunnel at Z = —107 m. 

3 Event reconstruction and selection 

A three-level trigger system was used to select events online [25,29]. At the first two levels, 
general characteristics of photoproduction collisions were required and background from 
beam-gas events was rejected. At the third level, jets were reconstructed by applying 
a cone jet algorithm to the CAL cells. Events with at least two jets satisfying E£ > 
4.5 GeV and if et < 2.5 were accepted, where E^ and if et are the transverse energy and 
pseudorapidity of the jet evaluated in the laboratory frame. 

In the offline analysis, the hadronic final state was reconstructed using energy-flow ob- 
jects [30,31] (EFOs), which are formed from a combination of track and calorimeter infor- 

2 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. 
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mation. This approach optimises the energy resolution and the one-to-one correspondence 
between the detector-level objects and the hadrons. The EFOs were corrected [32,33] to 
account for energy losses in the dead material and were forced to be massless by setting 
the energy component equal to the magnitude of the three-momentum. 

The jets were reconstructed using the kj< cluster algorithm [34] in the longitudinally 
invariant inclusive mode [35] using the p\ recombination scheme. At the detector level, 
jets were formed from the EFOs. The energy and angular reconstruction of the detector- 
level jets exhibited no systematic bias with respect to those defined at the hadron level. 

The data sample was selected by requiring the following: 

• the longitudinal position of the reconstructed vertex was in the range |Z vtx | < 40 cm; 

• 2/jb > 0.2, to reduce the contamination from beam-gas events, where y^ is the 
Jacquet-Blondel estimator [36] of the inelasticity, y, the fraction of the incoming 
electron momentum carried by the photon; 

• no scattered electron was observed in the CAL with y e < 0.85, to remove background 
due to DIS, where y e is the electron-method estimator of y; 

• Vjb < 0.85, which further removed contamination from DIS; 

• P™ iss /y/E^ < 2 GeV 1 / 2 , where P^ iss is the missing transverse momentum and 

is the total transverse energy. This removed any background from charged current, 
cosmic and halo-muon events; 

• at least three (or four, depending on the specific sample) jets were found with \rf ct \ < 
2.4 and E2? > 6 GeV. 

Furthermore, the three- and four-jet analyses were conducted in low- and high-mass re- 
gions given by 25 < M n j < 50 GeV and M n j > 50 GeV, respectively, where M n j is the 
invariant mass of the n-jet system. 

After all of the selection criteria had been applied, the low- and high-mass, three-jet data 
samples had 291646 and 38098 events, respectively, while the equivalent four-jet data 
samples had 31533 and 12525 events. 

Photoproduction events are characterised by the low virtuality, Q 2 , of the quasi- real 
exchanged photon, where Q 2 is the negative of the photon four-momentum squared. The 
kinematics are specified by y and x 7 , the fraction of the photon momentum taking part in 
the interaction. The variable be approximated using the observable x° bs , defined 

as 

o bs _ Ei=i E^eM-vf) , , 

7 " 2yE e ' [l) 

where the sum runs over the n-jets considered in the event and E e = 27.5 GeV denotes 

the energy of the incoming electron. 
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An angular variable used to analyse three-jet events is ip^ [37], which, as shown in Fig. [21 
is the angle in the three-jet centre-of-mass frame between the plane containing the highest 
energy jet and the beam, and the plane containing the three jets. The beam is described 
by the three- vector p beam = p p — p e , where p p and p e are the momenta of the proton and 
electron beams, respectively. The variable may be written as 

/ , \ (Pbcam X pg) ■ (p 4 X p 5 ) 

cos(^ 3 ) = -j n T' ( 2 ) 

IPbeam X p 3 ||P4 X p 5 1 

where it is conventional to number the three jets, 3, 4 and 5, in order of decreasing energy. 
The ips angle reflects the orientation of the lowest-energy jet. In the case where this jet 
arises from initial-state gluon radiation, the coherence property of QCD will tend to orient 
the third jet close to the incoming proton or photon direction. The two planes shown in 
Fig. [2] will therefore tend to coincide leading to a i(} 3 distribution that peaks toward and 

7T. 

4 Monte Carlo models 

Two MC generators were used to simulate photoproduction events, Herwig 6.505 [38] 
and Pythia 6.206 [39]. Both models include the LO (2 — > 2) matrix elements, approxi- 
mate higher-order processes using initial-state and final-state parton showers and simulate 
hadronisation. Direct and resolved photoproduction samples were generated separately. 
The implementation of the parton-showers and hadronisation models in the generators 
differs [39-42]. Hadronisation in Herwig is simulated using the cluster model [43] while 
Pythia uses the Lund string model [44,45]. 

4.1 Underlying-event models 

In the Herwig samples, MPIs were simulated using a separate program called Jimmy 
4.0 [46-49]. The latter is based on a simple eikonal model that approximates the inter- 
acting hadrons as disks that are extended in the transverse plane. An impact parameter 
quantifies the degree to which the two disks overlap during the collision. 

The Pythia MPI model that was used in this paper is known as the "simple model" [39]. 
It estimates the average number of MPIs per event, n, as n = a#(j3™ m )/crND(<s), where an 
and <tnd are the inclusive-hard and non- diffract ive, inelastic cross sections, respectively, 
and is a minimal constraint applied to the pt of the partonic collision. The number 
of MPIs is assumed to obey a Poisson distribution. The probability that the secondary 
collisions occur with transverse momentum is given by P(pr) = (V^ND^Xdcrtf/dpr)- 



4 



The simple model proceeds by ordering the secondary interactions in terms of pt and then 
calculates the probability of each successive scatter. 

4.2 Monte Carlo samples 

The MC samples were used for two purposes: to correct the data for detector effects and 
to compare to the measured hadron-level cross sections. The MC samples used to correct 
the data included a full Geant 3.13 [50] simulation of the ZEUS detector and three- level 
trigger. The resolved and direct MC samples were combined in the ratio that gave the 
best fit to the x° bs distribution in the data, significantly improving the overall description 
of the data by the MC. The resolved and direct samples used to compare to the measured 
cross sections, however, were combined in the ratio predicted by the MC models. 

The Herwig samples that were compared to the data were generated with ]3™ m set to 
4 GeV and with the CTEQ5L [51] and GRV-G LO [52] parameterisations for the proton 
and photon parton density functions (PDFs), respectively. All other Herwig parameters 
were set to default. To simulate the Herwig MPIs, the Jimmy model was run with the 
minimum p T of the secondary scatters, p™ pi , set to 2.2 GeV, the probability that the 
photon would resolve via a large hadron-like fluctuation set to 1/340 and the effective 
transverse radius of the photon squared was increased by a factor of 4 from its default 
value. These adjusted MPI settings were the result of simultaneous tuning to the three- 
and four-jet x° bs cross sections measured in this analysis, in the low- and high-mass 
regions. 

The Pythia samples that were compared to the data were generated with p™ m = 4.5 GeV, 
p^ pi = 1.9 GeV, which is the default value, and with the CTEQ5L and GRV-G LO 
parameterisations for the proton and photon PDFs, respectively. All other Pythia and 
simple-model parameters were set to default. 

4.3 Monte Carlo scaling factors 

Both LO MC simulations underestimated the magnitude of the measured cross sections. 
To compare with the shape of the differential cross sections, the MC predictions were 
either area normalised to the equivalent data set or scaled by the factor required to 
area normalise the high-mass MC cross sections to the high-mass data, where MPIs are 
expected to be less influential. The choice of normalisation procedure will be clearly 
stated in Section [8] and in the figure captions. 

The scaling factors that were applied to the Herwig and Pythia cross sections are shown 
in Table [TJ The scaling factors for the Pythia model were significantly larger than for 
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Herwig. Within each model, the factors applied to the three-jet predictions were similar 
whether MPIs were included or not, indicating that both models expect relatively little 
effect from MPIs in the high-mass three-jet sample. In the four-jet case, however, the 
scaling factors for the high-mass MPI samples were significantly smaller than those for 
the nominal samples, showing that both models predict a sizable contribution from MPIs 
in the high-mass four-jet sample. 



5 The 0{aa 2 s ) pQCD calculation 

The three-jet differential cross sections were calculated at O(ao 2 s ) using the program by 
Klasen, Kleinwort and Kramer [53,54]. The calculation is LO for this process. The 
renormalisation, /x r , and factorisation, fif, scales were set to // r = = = _E™ ax , the 
Et of the hardest parton. The theoretical scale uncertainty was evaluated by varying 
[A, setting it to 2 ±1 £'™ ax . The value of a s was calculated with one- loop precision and 
assuming five active flavours. A value of = 181 MeV was taken. The CTEQ6L [55] 
proton and GRV-G LO [56] photon PDFs were used for the calculation. 

The theoretical calculations were corrected for both hadronisation and MPI effects. The 
correction factors were obtained using the Herwig and Pythia models. The hadronisa- 
tion corrections, Chad, were calculated by taking the bin- by-bin ratio of the MC cross sec- 
tions at the hadron (ohl) an d parton levels (opl), Chad = chl/opl- To obtain cr PL , the kx 
cluster algorithm was run over all partons produced by the parton shower, prior to hadro- 
nisation. The MPI corrections, Cmpi, were calculated by taking the bin-by-bin ratio of the 
hadron-level cross sections, with (^ L PI ) and without MPIs K° L MPI ), C M pi = < l P 7^h°l MPI - 
The hadronisation and MPI corrections that were applied to the calculation represented 
the average corrections taken from the two MC models. A symmetric uncertainty was 
associated with the correction factors equal to half the difference between the predictions 
of the two MC models. The hadronisation corrections were found to typically reduce 
the cross section by 0(20%), whereas the MPI corrections increased the cross section by 
(9(50%) in the low-mass region and by (9(10%) in the high-mass sample. The uncertainty 
associated with the MPI corrections was larger than that estimated for the hadronisation 
factors. 



6 Acceptance corrections 

Detector effects and trigger inefficiencies were accounted for by applying bin-by-bin cor- 
rection factors, Cdet = 0"hl/0"dl> where cxdl is the predicted detector-level cross section. 
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The correction factors were extracted separately from the Herwig and Pythia MC sam- 
ples with MPIs. Before calculating the corrections, the resolved components of the MC 
samples were reweighted in to improve the overall description of the data [57]. The 
typical bin-by-bin corrections applied to the low- and high-mass three-(four-)jet data were 
approximately 0.75 (0.45) and 1.3 (0.95), respectively. The largest correction factors were 
found in the lowest and the extreme rf et bins. 

Once corrected to the hadron level, the measured cross sections, cr^g, at y/s = 300 GeV 
and y/s = 318 GeV were combined using the following formula: 

°300 " (°318/ a 30o) ' ^300 + °318 " £318 /0 \ 

° = 7 — T? ' (3) 

1-300 + ^318 

where C^yg is the integrated luminosity and er"^, the predicted cross section at y/s. As 
such, the measured cross sections presented here correspond to y/s = 318 GeV. The 
er^g/a^Q ratios were ~ 1.1 for all but the low-mass four-jet sample, where it was around 
10% larger. 

The results presented here represent the average hadron-level cross sections obtained 
when the data is treated with Herwig or Pythia, with half the spread interpreted as 
symmetric systematic uncertainty. This systematic uncertainty was added in quadrature 
to those described below. 



7 Systematic uncertainties 

A detailed study of the sources of systematic uncertainty associated with the measurement 
was performed using the Herwig MC sample [57]. The sources contributing to the quoted 
systematic uncertainties are, 

• the CAL energy scale uncertainty is ±3%. The energy scale was changed in the MC 
simulation accordingly; 

• the uncertainty due to the dependence of the measured cross sections on the detector- 
level selection criteria, which were varied up and down by the detector resolution for 
the data and the MC samples together. More specifically, 

— the Ej^ cut on each of the jets was changed by ±1 GeV; 

— the |?7j et | cut on each of the jets was changed by ±0.04; 

— the M n j selection criteria were changed by ±10%; 

— the lower and upper cut on y JB was changed by ±0.03 and ±0.05, respectively; 
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— the cut on y ej used to differentiate between real and fake scattered-electron can- 
didates, was changed by ±0.05. 

In addition, it was verified that the bin-by-bin acceptance corrections were not sensitive 
to variations in the relative amount of direct and resolved events or the exact shape of 
the i/jb distribution in the combined MC sample. 

All the systematic uncertainties were added in quadrature except that associated with 
the CAL energy scale. The CAL energy scale uncertainty is highly correlated between 
bins and is displayed separately in the plots presented here. 

The largest systematic uncertainty in the cross sections, except for the high-mass three- 
jet case, was that associated with the CAL energy scale, which led to an uncertainty 
of approximately ±10% in both the low- and high-mass three-jet cross sections, and 
±20% and ±15% in the low- and high-mass four-jet cross sections, respectively. The 
largest uncertainty in the high-mass three-jet cross section came from the choice of MC 
model used to calculate the acceptance corrections and was approximately ±10%. In the 
low-mass three-jet sample, the uncertainty due to the MC model was ±5% and in the 
four-jet low- and high-mass samples the values were ±3% and ±4%, respectively. Another 
significant source of systematic uncertainty in the low-mass samples was that associated 
with the Ej?" selection criteria, which generated a ±4% and ±8% effect in the three- and 
four-jet cases, respectively. In the high-mass samples, varying the M n j selection criteria 
generated a ±8% and ±6% effect in the three- and four-jet cross sections, respectively. 
The high-mass four-jet cross section was also sensitive to varying the E^ cut, which 
changed its value by about ±4%. 



8 Results and discussion 

The three- and four-jet photoproduction cross sections are presented here at the hadron 
level for jets with E^ > 6 GeV and |77-> et | < 2.4, in the kinematic region given by 
Q 2 < 1 GeV 2 and 0.2 < y < 0.85, and in low- (25 < M nj < 50 GeV) and high-mass 
(M nj > 50 GeV) regions. 



8.1 The da/dM n j cross section 

The three- and four-jet cross sections are given as a function of M n j in Fig. [3] and Tables [2] 
and [31 In general, both cross sections decrease exponentially with increasing M n j. Also 
shown in Fig. [3]are the Herwig and Pythia predictions without MPIs, normalised to the 
high-mass region (M nj > 50 GeV). Both models incorrectly describe the M n j dependence 
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of the cross section and significantly underestimate the low mass data. The discrepancy 
is larger in the four-jet case. With the inclusion of MPIs, both scaled MC predictions give 
a reasonably good description of the data over the full M n j ranges. 



8.2 The dcr/d#° cross section 

The three- and four-jet cross sections are given as a function of x° bs for the low- and high- 
mass samples in Fig. H] and Tables H] and The distributions exhibit a peak at x° bs ~ 0.9 
in all but the low-mass four-jet sample. The low-mass cross sections show enhancement at 
low x° bs compared to those at high-mass, which is a consequence of the tighter kinematic 
constraints at high mass. The degree to which the low x° bs region is enhanced is larger 
in the four-jet case. Shown also in Fig. @]are the scaled MC predictions compared to the 
data. 

The MC predictions show that, in the absence of MPIs, the three- and four-jet cross 
sections are expected to decrease with decreasing x° bs in both mass regions. Moreover, 
the low-x° bs suppression is expected to be marginally stronger in the four-jet case. The 
low-mass data, however, are in contradiction to both predictions. These data therefore 
suggest that some mechanism in addition to the processes modelled in the MC without 
MPIs, is contributing to the low-mass cross sections. One possible mechanism is MPIs. 

A comparison of the data with the MC predictions, normalised to the high-mass data, 
shows that both those with and without MPIs are in reasonable agreement with the high- 
mass cross sections (Figs. |4b andHU). In the low- mass low-x° bs regions that are poorly 
described by the nominal MC models, the introduction of MPIs to the simulations aids 
the description of the data. More specifically, the Herwig prediction with tuned MPIs 
describes the data reasonably well in all of the samples. The default Pythia model tends 
to overestimate the data in the mid-x° bs region (0.4 < x° bs < 0.8) in the low- mass samples 
but describes the cross section reasonably well elsewhere. 

Shown also in Fig. [4] is the contribution to the Herwig predictions attributed to direct 
processes, as defined in that LO model with parton showers. In all four samples, direct 
events are predicted to be found solely at high x° bs , although resolved processes are 
expected to make a significant, and in the four-jet case dominant, contribution even at 
high x° bs . This is in contrast to the behaviour of dijet photoproduction, where the high- 
x° bs region is dominated by direct events [15]. 



9 



8.3 The da/dy cross section 



The three- and four-jet cross sections are given as a function of y in the low- and high-mass 
samples in Fig. [5] and Tables [6] and [7J The cross sections are observed to increase steeply 
with y in all but the low-mass three-jet sample where it is approximately constant. The 
shapes of the distributions are governed by the available kinematic phase space and the 
energy distribution of the photon flux: while the phase space increases with y, the photon 
flux decreases with increasing y. 

Also shown in Fig. [5] are the predictions from Herwig and Pythia, with and without 
MPIs, which have been area normalised to the data. At high mass, the predictions with 
and without MPIs are similar and describe the data well. In the low-mass samples, MPIs 
are predicted to cause a more steeply increasing cross section. However, the data are not 
precise enough to differentiate between the various models and all of the MC predictions 
roughly describe the data. 

8.4 The da/dE^ and d a/drf ci cross sections 

The cross sections are given as a function of the of each of the jets in both three-jet 
samples in Fig. [6] and Tables [HI to [TU1 and in both four-jet samples in Fig. [7J and Tables [TT1 
to [141 Similarly, the da/drj }et cross sections are given in Figs. [HJ and M and Tables [151 to 
I2T1 In each case the jets have been ordered in descending ■ 

The E^ distributions of all of the jets in all of the samples are observed to fall off approx- 
imately exponentially with increasing E^ . All of the samples are similarly distributed 
when binned in terms of rj 3et , regardless of the position of the jet in the E 3 ^ ordering. The 
generic trend is an increasing cross section from rf ct « —1.4 followed by a plateau and, 
in the high-mass samples, a continued rise above rf et m 1.4. In the low- mass samples, 
the plateaux begin at rf et a; 0.6 and extend to the upper edge of the measured r\ region. 
In the high-mass samples, the lower bounds of the rf et ranges covered by the plateaux 
decrease with decreasing E^ , while the upper bounds all lie at rf ct 1.4. 

Shown also in Figs. [6] and [7J are da/dE^ predictions from the Pythia and Herwig 
models, with and without MPIs, normalised to the data. The description of the data 
by each of the MC models is generally good for all of the jets in all of the samples, 
although the MC predictions do vary somewhat. In all four samples, the PYTHIA cross 
sections without MPIs give the poorest description, in particular of da/dE 3 ^" 2 . Overall, 
the description of the data is improved by the inclusion of MPIs. 

Figures [HJ and [9] show the comparison between da/drf^ in the data and as predicted by 
the MC models. Each MC model largely describes both sets of high-mass data. The 
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best description of the low-mass data is given by the Herwig model with tuned MPIs. 
The Pythia prediction with default MPIs also generally describes the data. The poorest 
description is of the three-jet rf ctl distribution. The introduction of MPIs into the sim- 
ulations certainly aids the description. In general, the MC models without MPIs predict 
a cross section that falls off at high rf ct but MPIs reduce this effect in line with what is 
observed. 

8.5 The dcr/d cos(^3) cross section 

The three-jet cross section is given as a function of cos(t/>3) for the low- and high-mass 
samples in Fig. [TU] and Table [221 The cos(^3) cross section has a similar shape in both the 
low- and high-mass samples; relatively flat in the central cosf^) region and increasing 
rapidly as | cos^s)] — > 1. Also shown in Fig. [10] is a comparison of the MC models with 
the data. All of the area normalised MC predictions describe the data well. 

8.6 The da/dM^j cross section compared with 0{aa 2 s ) pQCD 

Figure [TT1 shows an 0{aa 2 s ) prediction, corrected for hadronisation effects and MPIs, com- 
pared to the measured da/dM 3 j cross section. The hadronisation and MPI corrections, 
including their estimated uncertainties, are given in Fig. [TTb and in Table [2j The hadro- 
nisation corrections are constant in M 3j -, while the MPI corrections increase significantly 
towards low M$j. The theoretical uncertainties on both the MPI corrections and the 
pQCD predictions are large. The magnitude and shape of the calculation is consistent 
with the data within the large theoretical uncertainties. This is best seen in the data over 
theory ratio shown in Fig. [TTb . The level of consistency between data and theory would 
be far worse at low M$j if it were not for the large MPI corrections. 

8.7 The da/dE^ cross section compared with 0(aa 2 s ) pQCD 

Figure [12] shows the comparison between the measured and predicted da/dE 3 ^ cross 
sections at 0(aa 2 s ) for each jet in the low- and high-mass three-jet samples. The pQCD 
predictions have been corrected for hadronisation effects and MPIs, which are detailed in 
Tables [51 to [TU1 (not shown in Fig. [12]). Hadronisation is predicted to decrease the parton- 
level cross sections at low E^ , whereas the MPIs are expected to do the reverse in the 
low-mass region: at high mass they are small. 

The calculation is largely consistent, within the large uncertainties, with both the low- 
and high-mass data. The exception to this is the rapid predicted falloff of the E 3 ^ 1 cross 
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section at high E^ 1 in the low-mass sample, which is not observed in the data. This 
effect is due to the three-parton kinematics within the theory, which precludes the highest 
momentum parton from carrying more than half of the total centre-of-mass energy. Thus, 
for the low-mass sample, the M n j < 50 GeV criterium translates into an Et < 25 GeV 
constraint for the partons in the theory, as indeed observed. 

8.8 The da / &rf ci cross section compared with 0{aa 2 s ) pQCD 

Figure [T3l shows the comparison between the measured and predicted da/drf et cross sec- 
tions at 0(aa 2 ) for each jet in the low- and high-mass three-jet samples. The hadroni- 
sation and MPI correction factors are given in Tables [15] to [T7] (not shown in Fig. [T3|) . 
Hadronisation is predicted to cause an overall decrease in the parton-level cross sections, 
most significantly at low rf et . The MPIs in the low- mass sample are predicted to cause a 
small increase in the cross section at low if et , becoming larger as if et increases. The cal- 
culation is consistent with the rf ct distributions. There is some indication of a difference 
in shape at lower masses, however, the uncertainties are again large. The description of 
the low-mass data would be worse if it were not for the MPI corrections. 

8.9 The dcr/d cos(^3) cross section compared with 0{aa 2 s ) pQCD 

The comparison between the measured and the predicted da/d costys) cross sections 
at 0(aa 2 ) is shown for the low-mass sample in Fig. O and for the high-mass region 
in Fig. [T5j The ratio of the low-mass data divided by the calculation is shown along 
with the hadronisation and MPI corrections, which are also given in Table [22j In both 
low- and high-mass regions, the QCD calculation is consistent with the data. There is 
some indication of a difference in shape at low masses; however, the theory has large 
uncertainties. 

9 Summary 

Three- and four-jet states have been measured in hard 7p collisions at HERA, using 
an integrated luminosity of 121 pb -1 . The three- and four-jet cross sections have been 
measured for jets with > 6 GeV and l^ *! < 2.4, in the kinematic region given by 
Q 2 < 1 GeV 2 and 0.2 < y < 0.85 and in two mass regions defined as 25 < M n j < 50 GeV 
and M n j > 50 GeV. The three-jet events have been measured with over seven times the 
luminosity of the previous ZEUS publication and in a wider Msj region. The four-jet 
process described here has been measured for the first time at HERA. 
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In the high-mass regions, the shape of the three- and four-jet cross sections are reasonably 
well described by both the Pythia and Herwig models without MPIs. In the low-mass 
region, the MC models without MPIs underestimate the data when normalised to the 
measured high-mass cross section. When MPIs are added to the MC simulations, the 
agreement between the models and data is generally improved. Although the data have 
large uncertainties, the measured cross sections are potentially useful in the testing and 
tuning of MPI and underlying-event models. 

The 0{aa 2 s ) pQCD calculation was compared to the three-jet data and is consistent within 
the large uncertainties. The MPI corrections typically improved the description of the 
data by the pQCD calculation. This data will provide a testing ground for higher-order 
calculations in photoproduction. 
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HERWIG scaling factors 


3-jet no MPIs 


3-jet with MPIs 


4-jet no MPIs 


4-jet with MPIs 


1.7 


1.7 


3.1 


2.1 


Pythia scaling factors 


3-jet no MPIs 


3-jet with MPIs 


4-jet no MPIs 


4-jet with MPIs 


3.8 


3.1 


9.2 


5.3 



Table 1: The scaling factors applied to the Herwig and Pythia x° hs and M, 
cross sections. 



M 3j 


range 


da/dM 3j 


±stat. 


±syst. 


±E-scale 




had. 


corr. 


MPI 


corr. 


(GeV) 






(pb/GeV) 
















25 


- 32 


94 


± 







+10 
-9 


+12 
-12 


0. 


793 


±0 


004 


1.88 


±0.37 


32 


- 41 


80.4 


± 





3 


+8.5 
-8.2 


+9.5 
-9.5 


0. 


805 


±0 


010 


1.44 


±0.23 


41 


- 50 


48.4 


± 





3 


+5.5 
-5.4 


+5.6 
-5.6 


0. 


815 


±0 


015 


1.25 


±0.16 


50 


- 59 


25.2 


± 





2 


+3.1 
-3.1 


+2.5 
-2.5 


0. 


815 


±0 


017 


1.15 


±0.13 


59 


- 70 


10.7 


± 





1 


+1.9 
-1.4 


+1.2 
-1.2 


0. 


817 


±0 


008 


1.09 


±0.10 


70 


- 82 


4.00 


± 





08 


+0.56 
-0.52 


+0.50 
-0.50 


0. 


810 


±0 


009 


1.07 


±0.08 


82 


- 95 


1.44 


± 





04 


+0.24 
-0.27 


+0.20 
-0.20 


0. 


808 


±0 


007 


1.04 


±0.05 


95 


- 109 


0.61 


± 





02 


+0.28 
-0.21 


+0.07 
-0.07 


0. 


811 


±0 


014 


1.03 


±0.05 


109 


- 124 


0.179 


± 





012 


+0.084 
-0.066 


+0.023 
-0.023 


0. 


840 


±0 


019 


1.03 


±0.05 


124 


- 140 


0.066 


± 





008 


+0.022 
-0.019 


+0.014 
-0.014 


0. 


798 


±0 


008 


1.00 


±0.04 


140 


- 160 


0.023 


± 





005 


+0.019 
-0.010 


+0.012 
-0.012 


0. 


809 


±0 


012 


1.00 


±0.06 



Table 2: The measured three-jet differential cross-section da/dM 3 j. The sta- 
tistical, systematic and calorimeter energy scale (E-scale) uncertainties are shown 
separately. Also shown are the hadronisation and MPI corrections that were applied 
to the 0{aa 2 s ) prediction. 
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M 4j 


range 


da/dM^j 


±stat. 


±syst. 


±E-scale 


(GeV) 






(pb/GeV) 




25 


- 32 


1.16 


± 


0.04 


+0.29 
-0.20 


+0.28 
-0.28 


32 


- 41 


5.8 


± 


0.1 


+0.6 
-0.7 


+1.1 
-1.1 


41 


- 50 


6.3 


± 


0.1 


+0.7 
-0.7 


+1.1 
-1.1 


50 


- 59 


4.49 


± 


0.09 


+0.65 
-0.48 


+0.73 
-0.73 


59 


- 70 


2.85 


± 


0.07 


+0.43 
-0.59 


+0.43 
-0.43 


70 


- 82 


1.25 


± 


0.05 


+0.27 
-0.20 


+0.22 
-0.22 


82 


- 95 


0.488 


± 


0.027 


+0.064 
-0.063 


+0.087 
-0.087 


95 


- 109 


0.237 


± 


0.018 


+0.052 
-0.068 


+0.042 
-0.042 


109 


- 124 


0.060 


± 


0.009 


+0.019 
-0.013 


+0.010 
-0.010 


124 


- 140 


0.027 


± 


0.006 


+0.019 
-0.015 


+0.006 
-0.006 


140 


- 160 


0.0024 


± 


0.0016 


+0.0053 
-0.0019 


+0.0012 
-0.0012 



Table 3: The measured four- jet differential cross-section da/dM^j. Other details 
as described in the caption to Table® 
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x° bs range 



da/dx. 



obs 

7 



±stat. ±syst. 
(Pb) 



±E-scale 







25 < M 3j 


< 50 GeV 


0.00 ■ 


- 0.11 


109 


± 7 


0.11 ■ 


- 0.22 


1449 


± 15 


0.22 ■ 


- 0.33 


2102 


± 15 


0.33 ■ 


- 0.47 


1989 


± 13 


0.47 ■ 


- 0.62 


1774 


± 12 


0.62 ■ 


- 0.75 


1752 


± 13 


0.75 ■ 


- 0.85 


2579 


± 19 


0.850 


- 0.935 


3803 


± 26 


0.935 


- 1.000 


730 


± 12 



+66 


+31 


-53 


-31 


+334 


+289 


-324 


-289 


+240 


+361 


-230 


-361 


+128 


+292 


-119 


-292 


+117 


+234 


-97 


-234 


+115 


+199 


-99 


-199 


+204 


+253 


-184 


-253 


+672 


+227 


-600 


-227 


+203 


+22 


-194 


-22 



M 3j > 50 GeV 



0.11 ■ 


- 0.22 


6.9 


± 1.0 


+8.0 
-5.0 


+2.0 
-2.0 


0.22 ■ 


■ 0.33 


51 


± 3 


+10 

-13 


+10 
-10 


0.33 ■ 


■ 0.47 


153 


± 4 


+42 
-42 


+21 
-21 


0.47 ■ 


■ 0.62 


310 


± 6 


+87 
-75 


+40 
-40 


0.62 ■ 


■ 0.75 


504 


± 8 


+93 
-86 


+60 
-60 


0.75 ■ 


■ 0.85 


871 


± 12 


+101 
-113 


+100 
-100 


0.850 ■ 


- 0.935 


1695 


± 18 


+202 
-211 


+161 
-161 


0.935 ■ 


■ 1.000 


666 


± 14 


+169 
-164 


+52 
-52 



Table 4: The measured three-jet differential cross-section da/dx^ 38 in both the 
low- and high-mass regions. Other details as described in the caption to Tabled 
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obs 


range 


da/dx 


° bs ±stat. 


±syst. 


±E-scale 








(Pb) 




25 < M Aj < 50 GeV 


00 

W . WW 


- 1 1 

W . ± ± 


4 9 


4- 1 7 


+5.1 
-3.0 


+1.4 
-1.4 


1 1 

W. ± J. 


- 99 


1 00 

X WW 


+ A 


+44 
-37 


+21 
-21 


0.22 


- 0.33 


210 


± 5 


+46 
-45 


+53 
-53 


0.33 


- 0.47 


177 


± 4 


+31 
-18 


+38 
-38 


0.47 


- 0.62 


136 


± 3 


+21 
-13 


+29 
-29 


0.62 


- 0.75 


94 


± 3 


+27 
-17 


+14 
-14 


0.75 


- 0.85 


113 


± 3 


+42 
-35 


+15 
-15 


0.850 


- 0.935 


125 


± 4 


+28 
-33 


+9 
-9 


0.935 


- 1.000 


22 


± 2 


+6 
-13 


+3 
-3 






M 4j 


> 50 GeV 






0.11 


- 0.22 


0.4 


± 0.2 


+2.3 
-0.2 


+0.2 
-0.2 


0.22 


- 0.33 


25 


± 2 


+14 
-12 


+9 
-9 


0.33 


- 0.47 


64 


± 3 


+26 
-23 


+11 
-11 


0.47 


- 0.62 


103 


± 3 


+26 
-30 


+19 
-19 


0.62 


- 0.75 


126 


± 4 


+18 
-24 


+21 
-21 


0.75 


- 0.85 


170 


± 5 


+31 
-27 


+24 
-24 


0.850 


- 0.935 


307 


± 8 


+39 
-40 


+35 
-35 


0.935 


- 1.000 


122 


± 5 


+69 
-37 


+7 
-7 



Table 5: The measured four-jet differential cross-section da/dx^ s in both the 
and high-mass regions. Other details as described in the caption to Tabled 
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y ran} 


?e 


da/dy 


±stat. 


±syst. 
(pb) 


±E-scale 




25 < M 3j 


< 50 GeV 




0.200 - 


0.272 


2264 


± 22 


+400 
-392 


+152 
-152 


0.272 - 


0.344 


2794 


± 22 


+491 
-382 


+286 
-286 


0.344 - 


0.417 


2854 


± 22 


+390 
-282 


+321 
-321 


0.417 - 


0.489 


2900 


± 22 


+288 
-238 


+340 
-340 


0.489 - 


0.561 


2973 


± 22 


+270 
-261 


+361 
-361 


0.561 - 


0.633 


2918 


± 22 


+268 
-242 


+363 
-363 


0.633 - 


0.706 


2833 


± 22 


+232 
-218 


+366 
-366 


0.706 - 


0.778 


2739 


± 22 


+250 
-216 


+380 
-380 


0.778 - 


0.850 


2797 


± 23 


+340 
-326 


+415 
-415 


M 3j > 50 GeV 


0.200 - 


0.272 


156 


± 7 


+24 
-38 


+16 
-16 


0.272 - 


0.344 


271 


± 7 


+28 
-13 


+24 
-24 


0.344 - 


0.417 


455 


± 10 


+78 
-65 


+43 
-43 


0.417 - 


0.489 


565 


± 11 


+59 
-91 


+52 
-52 


0.489 - 


0.561 


684 


± 12 


+134 
-108 


+67 
-67 


0.561 - 


0.633 


768 


± 13 


+ 111 
-131 


+81 

-81 


0.633 - 


0.706 


905 


± 14 


+171 
-199 


+100 
-100 


0.706 - 


0.778 


961 


± 15 


+201 
-206 


+113 
-113 


0.778 - 


0.850 


1058 


± 16 


+259 
-224 


+130 
-130 



Table 6: The measured three-jet differential cross-section da/dy in both the 
and high-mass regions. Other details as described in the caption to Tabled 
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y ran{ 


re 


da/dy 


±stat. 


±syst. 
(pb) 


±E-scale 




25 < M 4j 


< 50 GeV 




0.200 - 


0.272 


90 


± 4 


+30 
-25 


+8 
-8 


0.272 - 


0.344 


113 


± 4 


+23 
-14 


+18 
-18 


0.344 - 


0.417 


154 


± 5 


+28 
-24 


+29 
-29 


0.417 - 


0.489 


171 


± 5 


+30 
-17 


+33 
-33 


0.489 - 


0.561 


186 


± 5 


+15 
-23 


+37 
-37 


0.561 - 


0.633 


210 


± 5 


+30 
-27 


+39 
-39 


0.633 - 


0.706 


215 


± 5 


+31 
-28 


+42 
-42 


0.706 - 


0.778 


218 


± 6 


+23 
-24 


+44 
-44 


0.778 - 


0.850 


246 


± 7 


+37 
-43 


+52 
-52 


M Aj > 50 GeV 


0.200 - 


0.272 


16.8 


± 2.4 


+5.3 
-6.8 


+3.1 
-3.1 


0.272 - 


0.344 


40 


± 3 


+13 

-13 


+6 
-6 


0.344 - 


0.417 


81 


± 4 


+13 
-7 


+9 
-9 


0.417 - 


0.489 


108 


± 5 


+ 8 
-9 


+15 
-15 


0.489 - 


0.561 


136 


± 5 


+13 

-16 


+17 
-17 


0.561 - 


0.633 


171 


± 6 


+28 
-17 


+24 
-24 


0.633 - 


0.706 


228 


± 7 


+33 
-47 


+33 
-33 


0.706 - 


0.778 


266 


± 8 


+54 
-66 


+40 
-40 


0.778 - 


0.850 


266 


± 8 


+60 
-51 


+47 
-47 



Table 7: The measured four- jet differential cross-section da/dy in both the 
and high-mass regions. Other details as described in the caption to Tabled 
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zrijetl 

Ulrp 


range 


da/dEf a 


±stat. 


±syst 


±E-scale 




had. 


corr. 


MPI 


corr. 


(GeV) 






(pb/GeV) 


























25 < M 3j 


< 50 GeV 














6.0 


- 8.5 


107 


± 1 




+19 
-17 


+n 
-li 


0.744 


±0.026 


2.09 


±0.70 


8.5 


- 11.2 


248 


± 1 




+20 
-20 


+29 
-29 


0.791 


±0.012 


1.67 


±0.34 


11.2 


- 14.2 


169 


± 1 




+11 
-10 


+21 
-21 





819 


±0.009 


1.38 


±0.17 


14.2 


- 17.6 


75.2 


± 0.6 


+6.9 
-7.2 


+9.7 
-9.7 


0.842 


±0.003 


1.24 


±0.12 


17.6 


- 21.5 


19.9 


± 0.3 


+2.6 
-2.3 


+2.6 
-2.6 


0.834 


±0.017 


1.17 


±0.09 


21.5 


- 26.0 


2.54 


± 0.09 


+0.38 
-0.32 


+0.37 
-0.37 


0.774 


±0.064 


1.16 


±0.11 


26.0 


- 31.2 


0.200 


± 0.022 


+0.071 
-0.077 


+0.040 
-0.040 


0.725 


±0.080 


1.21 


±0.20 


M 3j > 50 GeV 


6.0 


- 8.5 


2.7 


± 





2 


+1.3 
-1.0 


+0.6 
-0.6 





515 


±0 


052 


1.42 


±0.29 


8.5 


- 11.2 


14.4 


± 





i 


+5.8 
-4.7 


+1.5 
-1.5 





646 


±0 


044 


1.28 


±0.20 


11.2 


- 14.2 


24.5 


± 





i 


+5.5 
-5.6 


+2.5 
-2.5 





759 


±0 


024 


1.17 


±0.16 


14.2 


- 17.6 


27.8 


± 





i 


+4.7 
-4.0 


+2.7 
-2.7 





839 


±0 


018 


1.13 


±0.13 


17.6 


- 21.5 


23.1 


± 





3 


+2.2 
-2.1 


+2.2 
-2.2 





878 


±0 


007 


1.09 


±0.09 


21.5 


- 26.0 


14.1 


± 





2 


+1.5 
-1.3 


+1.5 
-1.5 





864 


±0 


016 


1.07 


±0.08 


26.0 


- 31.2 


6.07 


± 





12 


+0.59 
-0.58 


+0.72 
-0.72 





859 


±0 


025 


1.06 


±0.07 


31.2 


- 37.0 


2.06 


± 





07 


+0.19 
-0.18 


+0.32 
-0.32 





848 


±0 


034 


1.04 


±0.05 


37.0 


- 43.9 


0.79 


± 





04 


+0.63 
-0.16 


+0.22 
-0.22 





838 


±0 


032 


1.02 


±0.04 


43.9 


- 51.8 


0.246 


± 





020 


+0.044 
-0.046 


+0.066 
-0.066 





847 


±0 


045 


1.04 


±0.07 


51.8 


- 60.8 


0.067 


± 





007 


+0.021 
-0.031 


+0.013 
-0.013 





881 


±0 


039 


1.00 


±0.04 


60.8 


- 71.0 


0.035 


± 





004 


+0.019 
-0.019 


+0.001 
-0.001 





870 


±0 


001 


1.02 


±0.07 



Table 8: The measured three-jet differential cross-section da/dE^ in both the 
low- and high-mass regions. Other details as described in the caption to Tabled 
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4 et2 


range 


dajdEf 2 


±stat. 


±syst 


±E-scale 




had. 


corr. 


MPI 


corr. 


(GeV) 






(pb/GeV) 


























25 < M 3j 


< 50 GeV 














6.0 


- 8.5 


383 


± 1 




+48 
-47 


+43 
-43 


0.772 


±0.023 


1.68 


±0.39 


8.5 


- 11.2 


229 


± 1 




+10 

-9 


+28 
-28 





829 


±0.010 


1.41 


±0.18 


11.2 


- 14.2 


59.0 


± 0.5 


+3.7 
-2.7 


+7.4 
-7.4 





868 


±0.018 


1.25 


±0.10 


14.2 


- 17.6 


10.4 


± 0.2 


+1.3 
-1.0 


+1.3 
-1.3 





869 


±0.005 


1.19 


±0.06 


17.6 


- 21.5 


0.84 


± 0.05 


+0.42 
-0.23 


+0.11 
-0.11 





861 


±0.042 


1.16 


±0.09 


M 3j > 50 GeV 


6.0 


- 8.5 


22.2 


± 





5 


+9.0 
-8.0 


+2.6 
-2.6 





626 


±0 


027 


1.26 


±0.22 


8.5 


- 11.2 


40 


± 


1 




+11 

-9 


+4 
-4 





759 


±0 


019 


1.16 


±0.15 


11.2 


- 14.2 


34.3 


± 





4 


+4.1 
-4.0 


+3.6 
-3.6 





861 


±0 


005 


1.11 


±0.12 


14.2 


- 17.6 


21.6 


± 





3 


+1.5 
-1.4 


+2.3 
-2.3 





891 


±0 


008 


1.08 


±0.09 


17.6 


- 21.5 


11.5 


± 





2 


+0.7 
-0.6 


+1.3 
-1.3 





898 


±0 


Oil 


1.07 


±0.07 


21.5 


- 26.0 


4.48 


± 





12 


+0.49 
-0.54 


+0.47 
-0.47 





888 


±0 


015 


1.06 


±0.07 


26.0 


- 31.2 


1.57 


± 





07 


+0.18 
-0.12 


+0.16 
-0.16 





881 


±0 


010 


1.06 


±0.07 


31.2 


- 37.0 


0.558 


± 





035 


+0.065 
-0.091 


+0.077 
-0.077 





892 


±0 


028 


1.03 


±0.05 


37.0 


- 43.9 


0.214 


± 





021 


+0.046 
-0.053 


+0.025 
-0.025 





898 


±0 


003 


1.01 


±0.03 


43.9 


- 51.8 


0.071 


± 





013 


+0.026 
-0.020 


+0.008 
-0.008 





898 


±0 


028 


1.00 


±0.05 


51.8 


- 60.8 


0.042 


± 





009 


+0.022 
-0.025 


+0.004 
-0.004 





765 


±0 


092 


1.02 


±0.10 



Table 9: The measured three-jet differential cross-section da/dE 1 ^ in both the 
low- and high-mass regions. Other details as described in the caption to Table® 
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E ] f t3 


range 


da/dE^ 3 


±stat. 


±syst. 


±E-scale 


had. 


corr. 


MPI 


corr. 


(GeV) 




(pb/GeV) 












25 < M 3j < 50 GeV 


6.0 


- 8.5 


644 


± 2 


+68 
-57 


+76 
-76 


0.799 


±0.011 


1.53 


±0.28 


8.5 


- 11.2 


63.2 


± 0.6 


+3.5 
-2.8 


+8.5 
-8.5 


0.859 


±0.011 


1.30 


±0.08 


11.2 


- 14.2 


3.82 


± 0.14 


+0.82 
-0.70 


+0.68 
-0.68 


0.887 


±0.017 


1.24 


±0.02 


14.2 


- 17.6 


0.096 


± 0.034 


+0.053 
-0.090 


+0.050 
-0.050 


0.834 


±0.058 


1.09 


±0.21 


M 3j > 50 GeV 


6.0 


- 8.5 


94 


± 1 


+20 
-20 


+11 

-11 


0.767 


±0.010 


1.15 


±0.14 


8.5 


- 11.2 


43.1 


± 0.5 


+5.1 
-1.1 


+4.4 
-4.4 


0.864 


±0.010 


1.08 


±0.09 


11.2 


- 14.2 


15.3 


± 0.3 


+1.5 

-1.3 


+1.7 
-1.7 


0.914 


±0.015 


1.06 


±0.07 


14.2 


- 17.6 


5.0 


± 0.1 


+1.0 
-1.0 


+0.5 
-0.5 


0.915 


±0.020 


1.04 


±0.05 


17.6 


- 21.5 


1.08 


± 0.06 


+0.15 
-0.15 


+0.16 
-0.16 


0.919 


±0.016 


1.04 


±0.06 


21.5 


- 26.0 


0.289 


± 0.026 


+0.086 
-0.068 


+0.034 
-0.034 


0.898 


±0.032 


1.00 


±0.04 


26.0 


- 31.2 


0.065 


± 0.012 


+0.030 
-0.034 


+0.006 
-0.006 


0.887 


±0.145 


1.00 


±0.08 


31.2 


- 40.0 


0.0185 


± 0.0053 


+0.0085 
-0.0078 


+0.0037 
-0.0037 


0.921 


±0.054 


1.05 


±0.15 



Table 10: The measured three-jet differential cross-section da/dE^ in both the 
low- and high-mass regions. Other details as described in the caption to Table [B 
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E^ 1 


range 


da/dE^ 1 


±stat. 


±syst. 


±E-scale 


(GeV) 




(pb/GeV) 








25 < M 4j 


< 50 GeV 






6.0 


- 8.5 


9.6 


± 0.2 


+1.4 
-1.6 


+1.7 
-1.7 


8.5 


- 11.2 


19.6 


± 0.3 


+1.9 
-2.0 


+3.7 
-3.7 


11.2 


- 14.2 


9.3 


± 0.2 


+1.4 
-1.1 


+1.8 
-1.8 


14.2 


- 17.6 


2.66 


± 0.10 


+0.40 
-0.52 


+0.50 
-0.50 


17.6 


- 21.5 


0.29 


± 0.03 


+0.12 
-0.06 


+0.06 
-0.06 


M 4j > 50 GeV 


6.0 


- 8.5 


1.00 


± 0.16 


+0.48 
-0.55 


+0.13 
-0.13 


8.5 


- 11.2 


6.1 


± 0.2 


+1.8 
-1.7 


+0.9 
-0.9 


11.2 


- 14.2 


8.3 


± 0.2 


+1.4 
-1.4 


+1.1 
-1.1 


14.2 


- 17.6 


6.46 


± 0.16 


+0.71 
-0.56 


+0.95 
-0.95 


17.6 


- 21.5 


3.93 


± 0.11 


+0.43 
-0.29 


+0.62 
-0.62 


21.5 


- 26.0 


1.93 


± 0.07 


+0.29 
-0.21 


+0.33 
-0.33 


26.0 


- 31.2 


0.74 


± 0.04 


+0.09 
-0.09 


+0.14 
-0.14 


31.2 


- 37.0 


0.289 


± 0.022 


+0.044 
-0.052 


+0.051 
-0.051 


37.0 


- 43.9 


0.100 


± 0.014 


+0.028 
-0.023 


+0.020 
-0.020 


43.9 


- 51.8 


0.016 


± 0.004 


+0.017 
-0.011 


+0.004 
-0.004 


51.8 


- 60.8 


0.0022 


± 0.0014 


+0.0040 
-0.0022 


+0.0010 
-0.0010 



Table 11: The measured four- jet differential cross-section da/dEj, in both the 
low- and high-mass regions. Other details as described in the caption to Tabled 



26 



rp range 
(GeV) 



da/dE^ t2 ±stat. ±syst. 

(pb/GeV) 



±E-scale 



6.0 

8.5 

11.2 

14.2 

17.6 

21.5 

26.0 

31.2 

37.0 



25 < M 4j < 50 GeV 



6.0 - 


8.5 


27.0 


± 0.3 


8.5 - 


11.2 


15.3 


± 0.3 


11.2 - 


14.2 


2.04 


± 0.10 


14.2 - 


17.6 


0.066 


± 0.017 



+3.1 

-3.6 

+2.2 

-1.8 

+0.50 

-0.30 

+0.050 

-0.061 



M 4j > 50 GeV 



8.5 

11.2 

14.2 

17.6 

21.5 

26.0 

31.2 

37.0 

43.9 



5.7 

12.7 

8.4 

3.40 

1.55 

0.55 

0.152 

0.058 

0.017 



± 0.2 
± 0.3 
± 0.2 
± 0.13 
± 0.08 
± 0.05 
± 0.021 
± 0.016 
± 0.005 



+1.9 
-1.5 
+1.8 
-1.4 

+0.5 

-0.7 

+0.48 

-0.25 

+0.26 

-0.51 

+0.17 

-0.09 

+0.043 

-0.032 

+0.036 
-0.029 
+0.014 
-0.009 



+4.9 

-4.9 

+2.9 

-2.9 

+0.33 

-0.33 

+0.004 

-0.004 



+0.7 
-0.7 
+1.9 
-1.9 
+1.4 
-1.4 

+0.48 

-0.48 

+0.24 

-0.24 

+0.08 

-0.08 

+0.036 

-0.036 

+0.013 

-0.013 

+0.005 

-0.005 



Table 12: The measured four- jet differential cross-section dajdE^ 2 in both the 
low- and high-mass regions. Other details as described in the caption to TablelM 
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E^ t3 range 
(GeV) 


da/dE^ ts ±stat. ±syst. ±E-scale 
(pb/GeV) 


25 < M 4j < 50 GeV 


6.0 - 8.5 
8.5 - 11.2 
11.2 - 14.2 


40.8 ±0.4 t¥ 7 

4.72 ±0.14 +_r 7 l ±8:1 
o.o54 ±0.019 +£88 +8;88| 


M 4j > 50 GeV 


6.0 - 8.5 
8.5 - 11.2 
11.2 - 14.2 
14.2 - 17.6 
17.6 - 21.5 
21.5 - 26.0 


17 c 1 n 4 +3.7 +2.6 
LI .O ZC U.^t _ 2 .6 -2.6 

13.1 ± 0.3 t\i t\Jj 

A a 1 Q 1 +0.7 +0.6 
31 u - i -1.5 -0.6 

n no ± 07 +0 - 21 +al ° 

31 u - u ' -0.22 -0.10 

0.191 ±0.031 ±8;8II 18:8£ 
0.035 ±0.009 tr ll ±8:881 



Table 13: TTie measured four- jet differential cross-section da/dE^ t3 in both the 
low- and high-mass regions. Other details as described in the caption to Tabled 



E^ u range 
(GeV) 


da/dE^ ±stat. ±syst. ±E-scale 
(pb/GeV) 


25 < M Aj < 50 GeV 


6.0 - 8.5 
8.5 - 11.2 


45.4 ± 0.4 tH ±8 6 

0.53 ±0.05 ±8if ±8:11 


M 4j > 50 GeV 


6.0 - 8.5 
8.5 - 11.2 
11.2 - 14.2 
14.2 - 20.0 


32.0 ±0.5 ttl ±49 
5.02 ±0.17 ±81 ±8£ 

o.8i ±0.07 ±8;^! ±8:88 
cose ±0.014 ±8;8^ 4 ±8:881 



Table 14: The measured four- jet differential cross-section da/dEj, in both the 
low- and high-mass regions. Other details as described in the caption to Tabled 
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range 


da/drf ctl 


±stat 


±syst. 


±E-scale 




had. 


corr. 


MPI 


corr. 










(pb) 
















25 < M 3j < 50 GeV 


-1.686 


- -1.371 


3.3 


± 


0.4 


+1.4 
-1.7 


+0.5 
-0.5 





324 


±0 


012 


1.15 


±0.21 


-1.371 


- -1.057 


36.0 


± 


1.4 


+7.8 
-9.9 


+2.9 
-2.9 





538 


±0 


041 


1.34 


±0.23 


-1.057 


- -0.743 


123 


± 


2 


+18 
-20 


+9 
-9 





686 


±0 


044 


1.33 


±0.18 


-0.743 


- -0.429 


256 


± 


3 


+40 
-37 


+19 
-19 





763 


±0 


025 


1.35 


±0.19 


-0.429 


- -0.114 


394 


± 


4 


+55 
-51 


+33 
-33 





793 


±0 


025 


1.35 


±0.19 


-0.114 


- 0.200 


506 


± 


4 


+67 
-63 


+51 
-51 





803 


±0 


014 


1.37 


±0.21 


0.200 
0.514 


- 0.514 

- 0.829 


600 

653 


± 
± 


5 
5 


+68 
-65 
+65 
-57 


+65 
-65 
+74 
-74 






820 
816 


±0 
±0 


009 
008 


1.39 
1.43 


±0.21 
±0.22 


0.829 
1.143 


- 1.143 

- 1.457 


631 
556 


± 
± 


5 
1 


+54 
-49 
+63 
-58 


+76 
-76 
+72 
-72 






814 
814 


±0 
±0 


004 
007 


1.49 
1.55 


±0.25 
±0.28 


1.457 


- 1.771 


573 


± 


4 


+58 
-50 


+85 
-85 





821 


±0 


003 


1.64 


±0.31 


1.771 


- 2.086 


697 


± 


6 


+97 
-84 


+103 
-103 





829 


±0 


002 


1.75 


±0.36 


2.086 


- 2.400 


790 


± 


9 


+186 
-186 


+107 
-107 





843 


±0 


005 


1.90 


±0.42 


M 3j > 50 GeV 


-1.686 


- -1.371 


2.4 


± 


0.5 


+1.7 
-1.6 


+0.6 
-0.6 





300 


±0 


040 


1.19 


±0.25 


-1.371 


- -1.057 


14.7 


± 


1.1 


+5.5 
-4.0 


+2.0 
-2.0 





498 


±0 


022 


1.20 


±0.15 


-1.057 


- -0.743 


42 


± 


2 


+12 
-11 


+4 
-4 





648 


±0 


030 


1.13 


±0.10 


-0.743 


- -0.429 


67 


± 


2 


+16 
-15 


+6 
-6 





757 


±0 


019 


1.10 


±0.10 


-0.429 


- -0.114 


84 


± 


2 


+18 
-16 


+7 
-7 





830 


±0 


019 


1.08 


±0.09 


-0.114 


- 0.200 


90 


± 


2 


+ 11 
-10 


+8 
-8 





863 


±0 


010 


1.06 


±0.06 


0.200 


- 0.514 


99 


± 


2 


+14 
-16 


+9 
-9 





867 


±0 


008 


1.06 


±0.07 


0.514 


- 0.829 


92 


± 


2 


+10 

-9 


+ 10 
-10 





840 


±0 


006 


1.07 


±0.08 


0.829 


- 1.143 


99 


± 


2 


+13 
-13 


+11 
-11 





830 


±0 


018 


1.08 


±0.09 


1.143 
1.457 


- 1.457 

- 1.771 


115 
145 


± 
± 


2 
2 


+14 
-19 
+17 
-21 


+14 
-14 
+18 
-18 






812 
808 


±0 
±0 


026 
019 


1.11 
1.12 


±0.12 
±0.12 


1.771 


- 2.086 


200 


± 


3 


+31 
-34 


+23 
-23 





818 


±0 


012 


1.16 


±0.16 


2.086 


- 2.400 


283 


± 


4 


+66 
-54 


+32 
-32 





829 


±0 


012 


1.22 


±0.17 



Table 15: The measured three-jet differential cross-section da/drf ctl in both the 
low- and high-mass regions. Other details as described in the caption to Tabled 
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^jet2 


range 


da/dr]i et2 


±stat 


±syst. 


±E-scale 




had. 


corr. 




MPI 


corr. 












(Pb) 


















25 < M 3j < 50 GeV 


-l.boo 


1 O 71 


y.i 


i 


O.O 


+7.4 
-3.4 


+1.8 
-1.8 





3oz 


_i_n 
±0 


005 


i 

1 


1 I 


_i_n on 
±U.Z(J 


-l.o/l 


1 n £7 


O / 


i 

± 


o 
Z 


+11 

-10 


+9 
-9 


U 


e>yz 


i n 
±U 


uiy 


i 

1 


6o 


i n o7 
±0.2 1 


-l.UO i 


n 7 A Q 


lOO 


± 


6 


+21 
-25 


+15 
-15 


n 
U 


1 zu 


±U 




i 

i 


i n 


±U.zo 


-0.743 


- -0.429 


282 


± 


1 


+39 
-30 


+26 
-26 





766 


±0 


015 


1 


42 


±0.26 


-0.429 


- -0.114 


408 


± 


1 


+41 
-44 


+39 
-39 





780 


±0 


012 


1 


43 


±0.25 


-0.114 


- 0.200 


492 


± 


1 


+54 
-46 


+50 
-50 





808 


±0 


007 


1 


42 


±0.22 


0.200 


- 0.514 


574 


± 


5 


+66 
-57 


+67 
-67 





817 


±0 


005 


1 


44 


±0.22 


0.514 


- 0.829 


628 


± 


5 


+67 
-61 


+73 
-73 





817 


±0 


003 


1 


45 


±0.21 


0.829 


- 1.143 


615 


± 


5 


+62 
-57 


+72 
-72 





816 


±0 


003 


1 


48 


±0.22 


1.143 


- 1.457 


550 


± 


4 


+40 
-34 


+66 
-66 





818 


±0 


013 


1 


50 


±0.25 


1.457 


- 1.771 


564 


± 


4 


+36 
-32 


+77 
-77 





819 


±0 


012 


1 


56 


±0.29 


1.771 


- 2.086 


686 


± 


6 


+79 
-82 


+96 
-96 





833 


±0 


015 


1 


63 


±0.33 


2.086 


- 2.400 


767 


± 


8 


+183 
-181 


+96 
-96 





844 


±0 


015 


1 


75 


±0.38 


M 3j > 50 GeV 


-2.000 


- -1.686 


0.48 


± 


16 +0 - 44 

u - iu -0.27 


+0.18 
-0.18 





085 


±0 


007 


1 


00 


±0.18 


-1.686 


- -1.371 


8.2 


± 


0.7 


+2.5 
-2.5 


+0.9 
-0.9 





347 


±0 


021 


1 


20 


±0.22 


-1.371 


- -1.057 


43 


± 


2 


+16 
-19 


+5 
-5 





571 


±0 


025 


1 


20 


±0.17 


-1.057 


- -0.743 


70 


± 


2 


+20 
-18 


+7 
-7 





710 


±0 


038 


1 


18 


±0.14 


-0.743 


- -0.429 


92 


± 


2 


+22 
-24 


+ 11 
-11 





783 


±0 


014 


1 


14 


±0.12 


-0.429 


- -0.114 


97 


± 


2 


+20 
-15 


+10 
-10 





822 


±0 


013 


1 


12 


±0.09 


-0.114 


- 0.200 


96 


± 


2 


+12 

-12 


+10 
-10 





845 


±0 


006 


1 


10 


±0.10 


0.200 


- 0.514 


104 


± 


2 


+12 
-13 


+11 
-11 





852 


±0 


007 


1 


09 


±0.08 


0.514 


- 0.829 


106 


± 


2 


+16 
-14 


+12 
-12 





858 


±0 


001 


1 


09 


±0.09 


0.829 


- 1.143 


108 


± 


2 


+14 

-17 


+ 11 
-11 





851 


±0 


009 


1 


09 


±0.10 


1.143 


- 1.457 


103 


± 


2 


+16 
-16 


+ 11 
-11 





846 


±0 


011 


1 


09 


±0.10 


1.457 


- 1.771 


113 


± 


2 


+19 
-14 


+13 
-13 





827 


±0 


018 


1 


10 


±0.11 


1.771 


- 2.086 


157 


± 


3 


+26 
-20 


+17 
-17 





829 


±0 


013 


1 


12 


±0.13 


2.086 


- 2.400 


232 


± 


4 


+51 
-49 


+23 
-23 





839 


±0 


010 


1 


14 


±0.15 



Table 16: The measured three-jet differential cross-section da/drf ct2 in both the 
low- and high-mass regions. Other details as described in the caption to Tabled 
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yyjet3 


range 


da/dr?** 3 


±stat. dbsyst. 
(pb) 


±E-scale 




had. 


corr. 


MPI 


corr. 


25 < M 3j < 50 GeV 


-2.000 
-1.686 


- -1.686 

- -1.371 


0.33 
15.2 


± 
± 


n 9 o +0.49 
u ' zo -0.33 

i.o Hi 


+0.08 
-0.08 
+1.3 
-1.3 






090 

388 


±0 
±0 


045 
019 


1 
1 


42 
30 


±0.82 
±0.33 


-1.371 


- -1.057 


80 


± 


o +16 
z -13 


+7 
-7 





631 


±0 


004 


1 


33 


±0.33 


-1.057 


- -0.743 


187 


± 


q +25 
° -24 


+17 
-17 





738 


±0 


007 


1 


38 


±0.34 


-0.743 


- -0.429 


308 


± 


4 +37 
^ -39 


+29 
-29 





791 


±0 


002 


1 


41 


±0.31 


-0.429 


- -0.114 


396 


± 


4 +46 
^ -35 


+42 
-42 





800 


±0 


006 


1 


44 


±0.30 


-0.114 


- 0.200 


474 


± 


4 +« 
1 _44 


+53 
-53 





826 


±0 


009 


1 


46 


±0.27 


0.200 
0.514 


- 0.514 

- 0.829 


560 
616 


± 
± 


c +67 
-60 

e +73 
O _74 


+67 
-67 
+72 
-72 






829 
818 


±0 
±0 


015 
000 


1 
1 


49 
51 


±0.26 
±0.22 


0.829 
1.143 


- 1.143 

- 1.457 


605 
538 


± 
± 


e +67 
-58 

4 +38 
^ -29 


+71 
-71 
+64 
-64 






807 
806 


±0 
±0 


010 
017 


1 
1 


52 
50 


±0.21 
±0.19 


1.457 


- 1.771 


559 


± 


4 +34 
^ -30 


+72 
-72 





808 


±0 


022 


1 


52 


±0.21 


1.771 


- 2.086 


685 


± 


c +79 
° -75 


+91 
-91 





822 


±0 


032 


1 


57 


±0.24 


2.086 


- 2.400 


793 


± 


o +183 
o _ 177 


+103 
-103 





835 


±0 


030 


1 


66 


±0.29 


M 3j > 50 GeV 


-2.000 


- -1.686 


1.9 


± 


0.3 t\l 


+1.0 
-1.0 





145 


±0 


014 


1 


01 


±0.11 


-1.686 


- -1.371 


26.1 


± 


i q +9.5 
L -° -6.9 


+3.1 
-3.1 





456 


±0 


004 


1 


19 


±0.22 


-1.371 


- -1.057 


70 


± 


o +14 
z -16 


+9 
-9 





673 


±0 


003 


1 


20 


±0.22 


-1.057 


- -0.743 


104 


± 


9 +26 
z -22 


+ 11 
-11 





788 


±0 


001 


1 


18 


±0.19 


-0.743 


- -0.429 


100 


± 


9 +20 
z -16 


+10 
-10 





840 


±0 


005 


1 


14 


±0.15 


-0.429 


- -0.114 


97 


± 


9 +12 
z -16 


+10 
-10 





876 


±0 


014 


1 


13 


±0.14 


-0.114 


- 0.200 


92 


± 


9 +11 
z -13 


+10 
-10 





901 


±0 


033 


1 


11 


±0.12 


0.200 


- 0.514 


89 


± 


9 +19 
z -15 


+9 
-9 





900 


±0 


036 


1 


10 


±0.11 


0.514 


- 0.829 


98 


± 


9 +17 
z -16 


+ 11 
-11 





856 


±0 


022 


1 


12 


±0.10 


0.829 
1.143 


- 1.143 

- 1.457 


95 
89 


± 
± 


9 +16 
z -14 

9 +11 

z -12 


+10 
-10 
+ 10 
-10 






840 
819 


±0 
±0 


002 
012 


1 
1 


11 
12 


±0.09 
±0.08 


1.457 


- 1.771 


103 


± 


9 +16 
z -15 


+12 
-12 





815 


±0 


023 


1 


11 


±0.06 


1.771 


- 2.086 


138 


± 


q +22 
° -21 


+16 
-16 





817 


±0 


030 


1 


11 


±0.08 


2.086 


- 2.400 


222 


± 


4 +44 
^ -49 


+19 
-19 





820 


±0 


042 


1 


12 


±0.07 



Table 17: The measured three- jet differential cross-section da/drf ct3 in both the 
low- and high-mass regions. Other details as described in the caption to Tabled 
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^jetl 


range 


da/dr]i etl 


±stat. 


±syst. 


±E-scale 








(pb) 








25 < M 4j 


< 50 GeV 






-1.371 


- -1.057 


0.61 


± 0.09 


+0.91 
-0.36 


+0.08 
-0.08 


-1.057 


- -0.743 


1.76 


± 0.23 


+0.85 
-0.67 


+0.30 
-0.30 


-0.743 


- -0.429 


6.1 


± 0.4 


+1.0 
-1.5 


+1.0 
-1.0 


-0.429 


- -0.114 


13.2 


± 0.6 


+2.2 
-3.0 


+2.1 
-2.1 


-0.114 


- 0.200 


21.9 


± 0.8 


+4.4 
-3.6 


+4.0 
-4.0 


0.200 


- 0.514 


36.7 


± 1.1 


+5.7 
-6.7 


+6.4 
-6.4 


0.514 


- 0.829 


43.3 


± 1.3 


+6.0 
-4.6 


+7.5 
-7.5 


0.829 


- 1.143 


55 


± 1 


+12 
-12 


+9 
-9 


1.143 


- 1.457 


43 


± 1 


+10 

-8 


+8 
-8 


1.457 


- 1.771 


44.5 


± 1.1 


+8.1 
-4.1 


+9.5 
-9.5 


1.771 


- 2.086 


54 


± 2 


+11 

-9 


+12 
-12 


2.086 


- 2.400 


68 


± 3 


+25 
-24 


+13 
-13 


M 4i > 50 GeV 


-1.371 


- -1.057 


1.10 


± 0.29 


+0.63 
-0.70 


+0.14 
-0.14 


-1.057 


- -0.743 


4.9 


± 0.6 


+2.3 
-1.4 


+0.6 
-0.6 


-0.743 


- -0.429 


10.9 


± 0.8 


+3.7 
-2.3 


+1.2 
-1.2 


-0.429 


- -0.114 


20.6 


± 0.9 


+5.1 
-7.1 


+2.5 
-2.5 


-0.114 


- 0.200 


20.9 


± 1.0 


+3.7 
-2.8 


+2.7 
-2.7 


0.200 


- 0.514 


24.0 


± 1.1 


+4.0 
-2.5 


+3.3 
-3.3 


0.514 


- 0.829 


24.4 


± 1.1 


+4.8 
-2.5 


+3.4 
-3.4 


0.829 


- 1.143 


27.8 


± 1.2 


+3.5 
-3.4 


+3.6 
-3.6 


1.143 


- 1.457 


28.0 


± 1.0 


+4.2 
-4.7 


+4.3 
-4.3 


1.457 


- 1.771 


37.0 


± 1.3 


+5.9 
-3.7 


+6.4 
-6.4 


1.771 


- 2.086 


47.0 


± 1.7 


+7.6 
-8.9 


+8.6 
-8.6 


2.086 


- 2.400 


63 


± 2 


+16 
-17 


+11 
-11 



Table 18: The measured four- jet differential cross-section da/drf ctl in both the 
low- and high-mass regions. Other details as described in the caption to Tabled 
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jet2 

'/ 


1 dll^c 


drrlrM ct2 

tiU I til I 


±stat. 


±syst. 


-\— n ,— Q C* Q 1 £i 








(pb) 










ou uev 






1 371 


-l.UO( 


n fi7 

U.O ( 


J. n 18 


+0.43 

-0.50 


+0.11 
-0.11 


1 0^7 


n 743 

-U. ( L iO 


0.10 


4- D V7 

ZLL U.O t 


+0.87 
-0.73 


+0.35 
-0.35 


n 743 

-U. ( ftO 


fl 49Q 


11.0 


± 0.5 


+6.5 
-6.9 


+1.0 

-1.0 


-D 429 


- -0 114 


1 7 


± 0.8 


+2.4 
-2.8 


+1.9 
-1.9 


-0.114 


- 0.200 


23.5 


± 0.9 


+4.1 
-3.6 


+3.5 
-3.5 


0.200 


- 0.514 


37.9 


± 1.2 


+4.9 
-5.8 


+6.3 
-6.3 


0.514 


- 0.829 


42.4 


± 1.3 


+4.9 
-5.7 


+7.5 
-7.5 


0.829 


- 1.143 


46.2 


± 1.3 


+4.2 
-6.4 


+9.0 
-9.0 


1.143 


- 1.457 


40.3 


± 1.1 


+8.7 
-6.9 


+8.3 
-8.3 


1.457 


- 1.771 


43.3 


± 1.1 


+8.7 
-5.5 


+9.0 
-9.0 


1.771 


- 2.086 


52 


± 1 


+ 10 

-7 


+11 
-11 


2.086 


- 2.400 


57 


± 2 


+19 
-15 


+11 
-11 


M 4i > 50 GeV 


-1.686 


- -1.371 


0.08 


± 0.04 


+0.21 
-0.08 


+0.03 
-0.03 


-1.371 


- -1.057 


3.1 


± 0.3 


+5.2 
-1.3 


+0.6 
-0.6 


-1.057 


- -0.743 


7.6 


± 0.7 


+1.9 
-1.4 


+1.0 
-1.0 


-0.743 


- -0.429 


15.2 


± 0.9 


+2.3 
-1.7 


+2.5 
-2.5 


-0.429 


- -0.114 


22.9 


± 1.2 


+4.2 
-6.0 


+3.5 
-3.5 


-0.114 


- 0.200 


20.9 


± 1.1 


+3.8 
-2.5 


+2.5 

-2.5 


0.200 


- 0.514 


26.4 


± 1.2 


+2.7 
-5.2 


+3.2 
-3.2 


0.514 


- 0.829 


33.8 


± 1.3 


+5.6 
-9.4 


+3.9 
-3.9 


0.829 


- 1.143 


28.1 


± 1.3 


+5.0 
-3.0 


+3.7 
-3.7 


1.143 


- 1.457 


25.4 


± 1.1 


+5.5 
-4.7 


+3.6 
-3.6 


1.457 


- 1.771 


29.4 


± 1.1 


+4.1 
-3.6 


+4.4 
-4.4 


1.771 


- 2.086 


43.3 


± 1.6 


+8.4 
-7.0 


+6.7 
-6.7 


2.086 


- 2.400 


51 


± 2 


+15 
-12 


+10 
-10 



Table 19: The measured four- jet differential cross-section da/drf ct2 in both the 
low- and high-mass regions. Other details as described in the caption to Tabled 
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dn/dr>i ct3 

LLU f til I 


±stat. 


±syst. 


-\— n Q C* Q 1 £i 








(Pb) 










\ JU vjrc V 






1 371 


1 n^7 




-i- n 1 

± U.l 


-0.4 


+0.3 
-0.3 


1 0^7 




3 7 


4- n a 


+1.1 

-0.9 


+0.9 
-0.9 


n 74^ 

-U. 1 fto 


ZL9Q 


Q A 


± 0.6 


+1.9 
-1.5 


+1.2 
-1.2 


-D 42Q 


- -0 114 


90 


± 0.8 


+7.2 
-7.8 


+3.2 
-3.2 


-0.114 


- 0.200 


26.2 


± 0.9 


+2.6 
-3.3 


+4.6 
-4.6 


0.200 


- 0.514 


37.2 


± 1.2 


+6.4 
-4.4 


+6.5 
-6.5 


0.514 


- 0.829 


45.7 


± 1.4 


+4.3 
-5.8 


+7.8 
-7.8 


0.829 


- 1.143 


47.8 


± 1.4 


+5.0 
-5.6 


+8.5 
-8.5 


1.143 


- 1.457 


37.8 


± 1.0 


+8.2 
-4.1 


+7.3 
-7.3 


1.457 


- 1.771 


40.0 


± 1.0 


+7.2 
-5.5 


+8.6 
-8.6 


1.771 


- 2.086 


51 


± 1 


+8 
-6 


+11 
-11 


2.086 


- 2.400 


61 


± 2 


+16 
-19 


+13 
-13 


M 4i > 50 GeV 


-1.686 


- -1.371 


0.59 


± 0.20 


+0.97 
-0.59 


+0.28 
-0.28 


-1.371 


- -1.057 


6.1 


± 0.6 


+2.0 
-2.0 


+1.4 
-1.4 


-1.057 


- -0.743 


19 


± 1 


+5 
-10 


+3 
-3 


-0.743 


- -0.429 


20.4 


± 1.1 


+4.6 
-4.3 


+4.6 
-4.6 


-0.429 


- -0.114 


21.0 


± 1.1 


+4.0 
-3.4 


+3.9 
-3.9 


-0.114 


- 0.200 


22.7 


± 1.2 


+4.8 
-4.3 


+3.6 
-3.6 


0.200 


- 0.514 


24.8 


± 1.2 


+6.0 
-4.1 


+3.7 
-3.7 


0.514 


- 0.829 


27.6 


± 1.3 


+3.8 
-3.4 


+4.1 
-4.1 


0.829 


- 1.143 


27.3 


± 1.3 


+6.3 
-4.1 


+4.1 
-4.1 


1.143 


- 1.457 


24.0 


± 1.1 


+4.1 
-2.7 


+3.9 
-3.9 


1.457 


- 1.771 


26.3 


± 1.0 


+4.3 
-2.7 


+3.9 
-3.9 


1.771 


- 2.086 


33.5 


± 1.4 


+4.3 
-2.9 


+4.8 
-4.8 


2.086 


- 2.400 


55 


± 2 


+11 
-13 


+5 
-5 



Table 20: The measured four-jet differential cross-section da/drf ct3 in both the 
low- and high-mass regions. Other details as described in the caption to Tabled 
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^jct4 


range 


da/drf eU 


±stat. 


±syst. 


±E-scale 








(pb) 




25 < M 4j < 50 GeV 


1 371 


-l.UO( 


i 3 


-1-0 9 


+1.0 

-i.i 


+0.3 
-0.3 


-1.057 


- -0.743 


4.4 


± 0.4 


+1.1 

-1.6 


+0.4 
-0.4 


-0.743 


- -0.429 


18 


± 1 


+10 
-13 


+3 
-3 


-0.429 


- -0.114 


18.8 


± 0.8 


+2.1 
-1.9 


+3.6 
-3.6 


-0.114 


- 0.200 


26.3 


± 0.9 


+4.0 
-2.4 


+4.5 
-4.5 


0.200 


- 0.514 


36.6 


± 1.2 


+4.5 
-4.3 


+6.4 
-6.4 


0.514 


- 0.829 


44.1 


± 1.3 


+7.1 
-6.3 


+7.7 
-7.7 


0.829 


- 1.143 


45.0 


± 1.4 


+7.3 
-5.0 


+7.9 
-7.9 


1.143 


- 1.457 


37.8 


± 1.0 


+7.9 
-4.1 


+7.0 
-7.0 


1.457 


- 1.771 


39.2 


± 1.0 


+7.5 
-3.5 


+8.2 
-8.2 


1.771 


- 2.086 


48.3 


± 1.3 


+7.6 
-6.8 


+9.8 
-9.8 


2.086 


- 2.400 


60 


± 2 


+14 
-15 


+12 
-12 


M 4i > 50 GeV 


-1.686 


- -1.371 


2.8 


± 0.4 


+2.0 
-1.8 


+5.6 
-2.2 


-1.371 


- -1.057 


7.3 


± 0.7 


+2.8 
-2.2 


+1.0 
-1.0 


-1.057 


- -0.743 


15.2 


± 0.9 


+3.8 
-2.2 


+2.0 
-2.0 


-0.743 


- -0.429 


21.6 


± 1.2 


+4.1 
-4.0 


+3.5 
-3.5 


-0.429 


- -0.114 


22.7 


± 1.2 


+5.5 
-5.7 


+3.5 
-3.5 


-0.114 


- 0.200 


24.3 


± 1.1 


+4.1 
-4.1 


+3.9 
-3.9 


0.200 


- 0.514 


23.3 


± 1.2 


+4.9 
-3.1 


+3.7 
-3.7 


0.514 


- 0.829 


27.3 


± 1.2 


+4.1 
-4.0 


+4.0 
-4.0 


0.829 


- 1.143 


34.6 


± 1.3 


+6.3 
-9.3 


+4.7 
-4.7 


1.143 


- 1.457 


25.8 


± 1.0 


+2.8 
-3.5 


+3.7 
-3.7 


1.457 


- 1.771 


25.6 


± 1.0 


+5.3 
-4.2 


+4.2 
-4.2 


1.771 


- 2.086 


34.0 


± 1.4 


+7.4 
-5.8 


+5.4 
-5.4 


2.086 


- 2.400 


44.3 


± 2.0 


+9.9 
-9.7 


+6.1 
-6.1 



Table 21: The measured four- jet differential cross-section da/drf ctA in both the 
low- and high-mass regions. Other details as described in the caption to Tabled 
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cos(ip 3 ) 


range 


da / dcos(tp3) 


±stat. 


±syst. 


±E-scale 




had. 


corr. 


MPI 


corr. 








(pb) 


















25 < M 3j < 50 GeV 


-1.00 - 
-0.88 - 


-0.88 
-0.72 


1494 
856 


± 13 
± 9 


+217 
-175 
+117 
-105 


+188 
-188 
+107 
-107 






864 
821 


±0.012 
±0.017 


1 
1 


76 
51 


±0.43 
±0.28 


-0.72 
-0.560 - 


-0.56 
-0.336 


804 
766 


± 8 
± 7 


+77 
-77 
+75 
-76 


+99 
-99 
+94 
-94 






793 
776 


±0.019 
±0.008 


1 
1 


43 
43 


±0.20 
±0.19 


-0.336 - 


-0.112 


759 


± 7 


+70 
-64 


+94 
-94 





763 


±0.005 


1 


42 


±0.20 


-0.112 - 


0.112 


746 


± 6 


+70 
-57 


+92 
-92 





761 


±0.006 


1 


41 


±0.19 


0.112 - 


0.336 


812 


± 7 


+80 
-80 


+98 
-98 





767 


±0.012 


1 


42 


±0.20 


0.336 - 
0.56 


0.560 
0.72 


851 
886 


± 7 
± 9 


+101 

-98 
+99 
-103 


+99 
-99 
+100 
-100 






780 
809 


±0.015 
±0.017 


1 
1 


48 
48 


±0.22 
±0.22 


0.72 


0.88 


898 


± 9 


+114 

-107 


+102 
-102 





839 


±0.013 


1 


51 


±0.27 


0.88 


1.00 


1475 


± 13 


+212 
-196 


+173 
-173 





894 


±0.006 


1 


64 


±0.38 










> 50 GeV 














-1.00 - 


-0.88 


542 


± 9 


+99 
-94 


+55 
-55 





824 


±0.018 


1 


15 


±0.15 


-0.88 - 


-0.72 


263 


± 5 


+44 
-45 


+29 
-29 





812 


±0.031 


1 


11 


±0.11 


-0.72 - 


-0.56 


215 


± 5 


+38 
-40 


+25 
-25 





807 


±0.021 


1 


10 


±0.10 


-0.560 - 


-0.336 


164 


± 3 


+26 
-24 


+19 
-19 





790 


±0.026 


1 


10 


±0.09 


-0.336 - 
-0.112 - 


-0.112 
0.112 


156 
149 


± 3 
± 3 


+31 
-30 
+33 
-27 


+18 
-18 
+17 
-17 






793 
785 


±0.012 
±0.013 


1 
1 


10 
11 


±0.10 
±0.11 


0.112 - 


0.336 


143 


± 3 


+25 
-22 


+15 
-15 





784 


±0.009 


1 


12 


±0.11 


0.336 - 


0.560 


151 


± 3 


+25 
-28 


+16 
-16 





802 


±0.017 


1 


12 


±0.11 


0.56 


0.72 


169 


± 4 


+32 
-34 


+17 
-17 





823 


±0.015 


1 


10 


±0.11 


0.72 


0.88 


217 


± 5 


+39 
-45 


+22 
-22 





847 


±0.006 


1 


10 


±0.11 


0.88 


1.00 


370 


± 7 


+67 
-61 


+38 
-38 





865 


±0.009 


1 


12 


±0.12 



Table 22: The measured three-jet differential cross-section da / dcos(ip 3 ) in both 
the low- and high-mass regions. Other details as described in the caption to Tabled 
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flow 



Figure 1: A schematic representation of an event exhibiting a secondary scatter 
(MPI). 



Pbeam 




plane (3,4,5) \ : 

---- plane (3,p bea J \\ 

Figure 2: A schematic representation of the three-jet centre- of-mass system and 
the 1JJ3 angle. 
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Figure 3: Measured cross section as a function of (a) M$j and (b) M 4j - for 
the three- and four- jet samples (solid circles). The inner error bars represent the 
statistical uncertainty. The outer error bars represent the statistical plus systematic 
uncertainties added in quadrature. The shaded band represents the calorimeter 
energy scale uncertainty. The predictions from the Herwig and Pythia models 
are also shown both with and without MPIs, as is the direct component of the 
Herwig cross sections. Each Monte Carlo cross section has been area normalised 
to the high-mass (M n ,j > 50 GeV) measured cross section by scaling the predictions 
by the factors indicated in the legend. 
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Figure 4: Measured cross section as a function of x° bs in the three-jet (a) low- 
and (b) high-mass samples and the four-jet, (c) low- and (d) high-mass samples. 
Other details as in the caption to Fig. 
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Figure 5: Measured cross section as a function of y in the three-jet (a) low- 
arid (b) high-mass samples and the four-jet, (c) low- and (d) high-mass samples. 
The predictions from the Herwig and Pythia models are also shown both with 
and without MPIs, as is the direct component of the Herwig cross sections. Each 
Monte Carlo cross section has been area normalised to the data. Other details as 
in the caption to Fig. 
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Figure 6: Measured cross section as a function of Eji in the three-jet (a) low- 
and (b) high-mass samples. The jets are ordered such that E^ f > E^ > E^ t3 . 
Each cross section has been multiplied by 10 n , where n is given on the plot, to aid 
visibility. Other details as in the caption to Fig. 
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Figure 7: Measured cross section as a function of E^ 1 in the four-jet (a) low- 
and (b) high-mass samples. The jets are ordered such that E^ 1 > E^ 1 > E^ t3 > 
E^ t4 . Each cross section has been multiplied by 10 n , where n is given on the plot, 
to aid visibility. Other details as in the caption to Fig. El 
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Figure 8: Measured cross section as a function ofrf et in the three-jet (a) low- and 
(b) high-mass samples. Each cross section has been shifted upwards by the amount 
given on the plot, to aid visibility. Other details as in the caption to Fig. 
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Figure 9: Measured cross section as a function ofif et in the jour-jet (a) low- 
and (b) high-mass samples. Other details as in the caption to Fig. 
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Figure 10: Measured cross section as a function o/cos(^3) in the three- jet (a) 
low- and (b) high-mass samples. Other details as in the caption to Fig. O 
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Figure 11: (a) Measured cross section as a function ofM 3 j in the three-jet sample 
compared with an 0(aa 2 ) prediction, corrected for hadronisation and MPI effects, 
(b) The hadronisation and MPI corrections as a function of M^j. (c) The ratio of 
the M 3 j cross section divided by the corrected 0{aa 2 s ) prediction. The theoretical 



uncertainty is represented by the dashed bands. 
FigM 



Other details as in the caption to 
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Figure 12: Measured cross section as a function of for each jet in the three- 
jet (a) low- and (b) high-mass samples compared with 0{aa 2 s ) predictions, corrected 
for both hadronisation and MPI effects. Each cross section has been multiplied by 
10 n , where n is given on the plot, to aid visibility. Other details as in the caption 
to Fig. [TO 
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Figure 13: Measured cross section as a function ofrf et for each jet in the three-jet 
(a) low- and (b) high-mass samples compared with 0{aa 2 s ) predictions, corrected 
for both hadronisation and MPI effects. Each cross section has been shifted upwards 
by the amount given on the plot, to aid visibility. Other details as in the caption 
to Fig. EH 
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Figure 14: (a) Measured cross section as a function of cos(^3) in the low-mass 
three- jet sample. Other details as in the caption to Fig. Q7J 
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Figure 15: (a) Measured cross section as a function o/cos(^>3) in the high-mass 
three- jet sample. Other details as in the caption to Fig. Q7J 
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